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ABSTRACT 

For a one-dimensional chain of four nuclear spins (1/2) and taking into account first 
and second neighbor interactions among the spin system, we make the numerical 
simulation of Shor prime factorization algorithm of the integer number A'^ = 4 
to study the influence of the second neighbor interaction on the performance of 
this algorithm. It is shown that the optimum Rabi's frequency to control the non- 
resonant effects is dominated by the second neighbor interaction coupling parameter 
(J'), and that a good Shor quantum factorization is achieved for a ratio of second 
to first coupling constant of J' / J > 0.04. 
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1. Introduction 



The polynomial time solution of the prime decomposition of an integer number, given 
by Shor factorization algorithm [1] using quantum computation, has triggered the 
huge amount of work in quantum computers [2] and quantum information [3] areas in 
physics. This algorithm has already been demonstrated for few qubits [4] quantum 
computers. A qubit is the superposition of two quantum states of the system, say 
|0) and |1), ^' = Co|0)+Ci|l) such that |CoP + |Cip = 1, and it is the basic element 
to process the information in a quantum computer. The states |0) and |1) can also 
be called basic qubits. The L-tensorial product of L-basic-qubits forms a register of 
length L, say \x) = \iL-i, ■■■ ,io) with ij = 0, 1 ("0" for the ground state and "1" for 
the exited state). The set of these states makes up the basis of the 2^-dimensional 
Hilbert space where the quantum computer works, and a typical element of this 
space is given by = J2 Cx\x) with X] |Ca;p = 1. A solid state quantum computer 
of our particular interest and which might be developed in a near future is using 
a one-dimensional chain of nuclear spins (1/2) which is inside a strong magnetic 
field in the z-direction and an rf-field in the x-y plane. The magnetic field in the z- 
direction determines the state of the nuclear spin, |0) if the nuclear spin is parallel to 
this magnetic field and |1) if the nuclear spin is in opposite direction. This magnetic 
field also determines the Zeeman spectrum of the system. The rf-ficld is used to 
cause the desired transitions between the Zeeman levels of the systems. Up to now, 
this model has been developed just theoretically and hopefully the technological and 
experimental part may start in a near future. However, because the Hamiltonian 
of this system is well known, very important theoretical studies have been made [5] 
which are also important for the general understanding of quantum computation 
[6]. In this model, first neighbor interaction among the nuclear spins is considered, 
and Shor quantum factorization of the number N = 4 has been simulated with 
this model [7]. In this paper, we consider also second neighbor interaction among 
the nuclear spins. Using this interaction, we study Shor factorization algorithm 
to factorize the integer number N=4, developing the proper code to do this. We 
study the well performance of this factorization through the fidelity parameter and 
determine the minimum value of the second neighbor interaction coupling constant 
to do this. Finally, we also point out the modification caused by the second neighbor 
interaction to the so called 27r/c-method, used to suppress non-resonant transitions. 

2. Equation of Motion 

Consider a one-dimensional chain of four equally spaced nuclear-spins system (spin 
one half) making an angle cos 6 = l/\/3 with respect the z-component of the mag- 
netic field (chosen in this way to kill the dipole-dipole interaction between spins) 
and having an rf-magnetic field in the transversal plane. The magnetic field is given 

by 

^ = (bcos{ujt + ip),—bsm{uot + ip),B{z)) , (1) 

where b, uj and (p arc the amplitude, the angular frequency and the phase of the 
rf-field, which could be different for different pulses. B{z) is the amplitude of the z- 
component of the magnetic field. Thus, the Hamiltonian of the system up to second 
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neighbor interaction is given by 

= - ^ • Bk - 2 ^ lUUi - 2^'^ E IkIU2 , (2) 

k=0 k=0 k=Q 

where /Xk represents the magnetic moment of the kth-nucleus which is given in terms 
of the nuclear spin as /Xk = h'y{I^,lf,I^), being 7 the proton gyromagnetic ratio. 
Bk represents the magnetic field at the location of the kth-spin. The second term 
at the right side of (2) represents the first neighbor spin interaction, and the third 
term represents the second neighbor spin interaction. J and J' are the coupling 
constants for these interactions. This Hamiltonian can be written in the following 
way 

H = Ho + W , (3o) 

where Hq and W are given by 

ifo = I E cufc/l + 2 J(7o^7f + If if + /|J|) + 2J'mil + nii) \ (36) 
U=o J 



and 

fiO 3 r 



2 ^ 

fc=0^ 



(3c) 



The term ojk represents the Larmore frequency of the kth-spin, = ^B{zk)- The 
term Q is the Rabi's frequency, J7 = 76. Finally, the term = ± ilf represents 
the ascend operator (+) or the descend operator (-). The Hamiltonian Hq is diagonal 
on the basis {|i3i2^i^o)}) where ij = 0, 1 (zero for the ground state and one for the 
exited state), 

-ffoK3«2n«o) = -Ej3i2jiioK3^2«l«o) • (4a) 
The eigenvalues Ei^i^i^i^ are given by 

Eisi^nio = -|{E(-l)^'=a;, + J ^(-l)^'=+^^+i + j' ^(-l)^'=+''=+A . (46) 

^fc=0 A;=0 k=0 ^ 



The term (3c) of the Hamiltonian allows to have a single spin transitions on the 
above eigenstates by choosing the proper resonant frequency. 

To solve the Schrodinger equation 

zh—=H^, (5) 
let us propose a solution of the form 

m = T.Ckm) , (6) 

where we have used decimal notation for the eigenstates in (4a), HQ\k) = E]^\k). 
Substituting (6) in (5), multiplying for the bra (m|, and using the orthogonality 
relation {m\k) = 6mki we get the following equation for the coefficients 

15 

ifiCm = EmCm + Y.Ck{m\W\k) m = 0, . . . , 15. (7) 

fe=0 
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Now, using the following transformation 

C„, = D^e-'""^''^ , (8) 

the fast oscillation term EmCm of Eq. (7) is removed (this is equivalent to go to the 
interaction representation), and the following equation is gotten for the coefficients 

Dm 

1 15 

iDm = TY. WmkDkc'^'-''' , (9a) 

"'k=Q 

where Wmk denotes the matrix elements (m|VF|A;), and Umk are defined as 

'^mk = ^ • (90) 

Eq. (9a) represents a set of 32 real coupling ordinary differential equations which 
can be solved numerically, and where W^k are given by 

Wmk = X (0, z "or" z*), (9c) 

where z is defined as z = e*('^*+'^) , and z* is its complex conjugated. 



3. Shor's factorization algorithm and numerical simulation 

Following Shor's approach to get the factorization of an integer number N, one 
selects a L + M-register of the form \x;y), where \x) is the input register of length 
L, and \y) is the valuation register of length AI. The \y) register will store the values 
of the periodic function y{x) = q^{mod N), where the integer g is a coprime number 
of N, that is, their grater common divisor is one (gcd{q,N) = 1). Thus, starting 
with the ground state, 

*o = |0;0), (10) 
the uniform superposition state is created in the x-register, 

*i = ^El^;o)' (11) 

X 

where 2^ is the number of qubits in the x-register. On the next step, the computation 
of the function y{x) = q^{mod N) is carried out in the y-register, 

X 

Then, the discrete Fourier transformation is done in the x-register, 

*3 = ^EE«'''^'^/''l^;2^(^))- (13) 

X k 

After these steps, one makes the measurement of the state on the x-rcgistcr, and 
the probability must be a peak distribution with peaks separation, Ax, equal to 
Ax = 2^/T, whenever 2^ be divisible by the period T. In this way, one finds the 
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period T of the function y{x) = q^{mod N). If this period is an even number, the 
factors of N can be computed finding the greatest common divisor of q^^'^ it 1 and 
the number N {gcd{q^/'^ ±l,N)). 

Now, for A^" = 4, one just needs two-qubits registers in the x-register (L = 2) and 
two-qubits registers in the y-register (M = 2). The only coprime number is g = 3, 
and the function y{x) = 3^ {mod 4) has period T = 2. Therefore, starting with the 
function 

$0 = |00;00) , (14a) 
the superposition state is created in the x-register, 

^r^ = 1||00;00) + |01;00) + |10;00) + |11;00)| . (146) 

Next, the function y{x) = 3^ {mod 4) is valuated on the y-register, 

*2 = ^||00;01) + |01;11) + |10;01) + |11;11)| . (14c) 

Then, the discrete Fourier transformation is performed to get (after summation of 
aU terms) the wave function 

*3 = ^||00;01) + |00;11) + |10;01) + 10;11)| . {14d) 

The measurement on the x-rcgister give us the states |00) or 10) (a; = or x = 2). 
So, one has Ax = 2, and the period of our function is T = 2^/2 = 2. Finally, the 
factors of iV = 4 (4 = 2 • 2) are obtained from gcd{3^/'^ -1,4) = 2. 

To make the numerical simulation of this algorithm, we have chosen the following 
parameters in units of 27rxMhz, 

ojQ = 100 , ui = 200 , U2 = 400 , W3 = 800, J = 10 , J' = 0.4 , = 0.1 . (15) 

These parameters were chosen in this way to have a clear separation on the Zeeman 
spectrum and to have a good definition for the resonant transitions in our numerical 
simulation. This does not imply a restriction on our simulations since our main 
results are applicable also to the current design [6]. Now, starting with the ground 
state of the system, |0000), we create a superposition state in the x-register using 
three 7r/2-pulses with zero phases and with resonant frequencies ljq^4, (Jq.s and 1^4,12- 
The valuation of the function y{x) = 3^ {mod 4) in the y-register is carried out with 
four TT-pulses with zero phases and with the resonant frequencies a;o,i, "^4,5, (^5,7 
and Wi3_i5. Finally, the discrete Fourier transformation in the x-register is gotten 
through five vr-pulses with zero phases and with the frequencies ujqj, (^2,6, ^2,3, 
'^14,15, and Willis- The transitions involved in the algorithm are shown in Fig. 1. 
Note that at the end of Shor's algorithm, one must get the following probabilities 

= ICsp = ICgp = |Cii|2 = 1/4 , (16) 

according with our final wave function (14d), note from (8) that |Cfc| = \Dk\. Fig. 2 
shows the behavior of the probabilities |Cfcp during the entire Shor's algorithm. Fig. 
2a shows the formation of the superposition state (14b), starting from the ground 
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state (14a). Fig. 2b shows the valuation of the function y{x) = 3^{mod 4) at the 
end of four vr-pulses, wave function (14c). Fig. 2c shows the formation of the wave 
function (14d) at the end of five 7r-pulses. To better illustrate what is happening 
during the Shor's algorithm, we calculated the expected values of the z-component 
of the spin of the system for each qubit. These expected values are given by 

m) = li:i-^f\Ck{t)\\ (17a) 

k=0 



{ID = ^||CoP + |Cip-|C2|2-|C3p + |C4p + |C5p-|C6|2-|C7p 



(176) 



{ID = ^||Co|' + |Ci|2 + |C2|2 + |C3|2-|C4|^-|C5|^-|C6|2-|C77p 

+|C8p + ICgp + ICiop + |Cnp - |Ci2p - |Ci3p - |Ci4p - ICisI 



(17c) 



and 



{ID = I j: \Ckf -if: \Ck\' . im 

k=0 k=8 

Fig. 3a shows these expected values during the formation of the superposition state 
on the x-register. Fig. 3b shows their behavior during the valuation of the function 
y{x) = 3^(mod 4) on the y-register, and Fig. 3c shows their behavior during the 
creation of the discrete Fourier transformation. Of course, this observed behavior is 
the behavior that one could have expected from the Shor's algorithm. Fig. 4a shows 
the probabilities of the expected four-qubits registers at the end of Shor's algorithm 
(wave function (14d)), and Fig. 4b shows the probabilities of the non-resonant 
states. 

To determine the value of the second neighbor coupling constant needed to have a 
good reproduction of the Shor's algorithm, we calculate the fidelity parameter [8], 

F = {^expected\^) , (18) 

where "^expected is the wave function (14d), and ^' is the resulting wave function 
from our computer simulation. This is done as a function of the ratio J'/ J (second 
to first neighbor coupling constant interactions). Fig. 5 shows our results of the 
calculation of this parameter as a function of J' /J. As one can see, for a value of 
J' / J > 0.04, one can say that we have already a very good reproduction of the 
Shor factorization algorithm. Of course, this does not mean that second neighbor 
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interaction is needed to implement Shor qiiantum algorithm sinee, if J' = 0, one can 
just change the protocol (resonant pulses) to get this algorithm. Rather, what this 
result means is that, in the case of having second neighbor interaction (J' 7^ 0), the 
optimum performance of our Shor quantum algorithm is gotten for the above ratio. 



4. Second neighbor interaction and the 27rfc-method 



One of the important results from the consideration of first neighbor interaction 
and the selection of the parameter as J/Ao; ^ 1 and fi/Aw <C 1 is the possibility 
of choosing the Rabi's frequency Vl in such a way that the non-resonant effects are 
eliminated. This procedure is called the 27rA;-method [9], and this Rabi's frequency 
for a TT-pulse is chosen as Q = |A|/\/4fc^ — 1, where k is an integer number, and A is 
the detuning parameter between the states \p) and |m), A = {Ep — Ejn)/h — uj being 
iu the electromagnetic resonant frequency. This detuning parameter is proportional 
to the first neighbor coupling constant J. Let us see how this detuning parameter is 
modified due to second neighbor interaction. Assuming that the states \p) and \m) 
are are the only ones involved in the dynamics, from Eq. (9a), one has 



it, 



and 



iD„ 



h 



(19) 



Thus, given the initial conditions -Dp(O) = Cp(0) and -Dm(O) = C'm(O)) the solution 
is readily given by 



Dp{t) = Cp{0) 



cos 



A nj, 

i— sin 

n, 2 



+ t — sm — - >e 2 (20a) 



and 



Dmit) = { Cm{0) 



cos ■ 



A Q.et 



JlCp(O) 
+ I ' sm 



n4 



■e 2 



(206) 



where Jig is defined as = VJ^^ + A^. So, for example, at the end of a 7r-pulse 
{t = T = 7r/Jl), the argument of the periodic functions in (20a) and (20b) is given 
by J7eVr/2J7. If one makes this term to be equal to any multiple of tt, one can get rid 
of the non-resonant terms since the solutions are given by 

Dp{T) = (-l)*^Cp(0)e^^-/2^ , and L>^(r) = (-l)^C^(0)e-^^-/2^ . 

The Rabi's frequency obtained from the condition $7e7r/2n = is given by J7 = 
|A|/\/4A;^ — 1, and this is the so called 27rA;-method. 

Now, let us select a resonant transition containing the Larmore frequency ujq. These 
frequencies could he cuq + J + J' , ujq — J + J' , lvq — J — J' or cuq + J — J' which 
correspond to the transitions (decimal notation) |0) |1) (|10) ^ |11)), |2) |3), 
1 6) 1 7), and |2) |3). So, all of these states are pertubated during a pulse, and the 
frequency difference A may have the values 2 J, 2 J', 2J + 2J', or 2J — 2J'. Choosing 
other Larmore frequencies, the additional values for the detuning parameter are 4J, 

(k) 

and 4J-I-2J . Thus, denoting by 0\ the Rabi's frequency selected by this method. 



n 



|A| 
V4fc^ - 1 



(21) 
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there are five possible values for A which are 4J + 2J', 4J, 2J + 2J', 2 J and 2 J'. 



To see the dependence of the Shor's algorithm with respect the Rabi's frequency, 
we use again the fidelity parameter parameter (18). Using the same values for our 
parameters as (15) but il, Fig. 6 shows the fidelity parameter as a function of the 
Rabi's frequency. The lines LI, L2 and L3 mark the fl's values where this fidelity 
has peaks. These peaks correspond to the following 27rA;-method Rabi frequencies 

^4J+2J' ^ ^4J ~ ^2,7+2,7' ^ ^27 ^ ^ ^27' ~ 0.080403 , 
^4J+2J' ~ ^4J^^ ~ ^2J+2J' ~ ~ ^2 J' = 0.100791 , 



and 



'47+2 J' ~ ^'4J ~ ''''2J+2J' ~ '''2J ~ "2 J' 

As one can see in Fig. 7, where we have plotted Jl^^ (for the detuning values 
mentioned above) and where the lines LI, L2 and L3 have been drawn, around 

these lines there are several other values of ^^4j'_|_27'' ^47^ ^2J+2J' ^^'^ ^2^J which, 
in principle, should cause a peak in the fidelity parameter (because they belong 
to the 27rA;-method) . However, they do not appear at all on Fig. 6. This means 
that the peaks values on the fidelity parameter are fully dominated by the second 
neighbor coupling interaction parameter (J'). 



5. Conclusions 

For a one-dimensional chain of nuclear spins (one half) system quantum computer, 
we have considered first and second neighbor interactions, and we have study the 
effects of second neighbor interaction on Shor quantum factorization algorithm of 
the number = 4. We have shown that a good factorization algorithm can be 
gotten if the ratio of second to first neighbor interaction constants is chosen such 
that J' /J > 0.04. We have shown also that the application of 27rA;-method to 
eliminate non-resonant transitions is not so simple since the detuning factor varies 
with both parameters J and (first and second neighbor coupling interactions). 
However, the peaks on the fidelity parameter are dominated by the second neighbor 
coupling parameter. In other words, if there exists second neighbor interaction in a 
chain of nuclear spin system quantum computer, this interaction may dominate the 
27rA;-method to control non-resonant effects. 
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Figure Captions 

Fig. 1 Energy levels and resonant frequencies used within the algorithm. 

Fig. 2 Probabilities |Cfcp : (k). (a) Formation of superposition state in the x- 
register, wave function (14b). (b) Formation of the wave function (14c). (c) Forma- 
tion the wave function (14d). 

Fig. 3 Expected values (/|): (k=0,l,2,3). (a) During formation of wave function 
(14b). (b) During formation of the wave function (14c). (c) During formation of 

the wave function (14d). 

Fig. 4 Probabilities |(7fcp. (a) For the expected states k = 1,3,5,7. (b) For the 
non-resonant states k = 0,2, 4, 6, 8, 9, 10, 11, 12, 13, 14, 15. 

Fig. 5 Fidelity parameter function of J' /J. 

Fig. 6 Fidelity parameter function of il. 

Fig. 7 Rabi frequency 17^^ as a function of A; for A = 4 J + 2J' (1), A = 4J (2), 
A = 2J + 2J' (3), A = 2J (4), A = 2J' (5). Lines LI, L2 and L3 correspond to 
Fig. 6. 
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